Universality class of quantum criticality for strongly repulsive spin-1 bosons with 

antiferromagnetic spin-exchange interaction 
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Using the thermodynamic Bethe ansatz equations we study the quantum phase diagram, thermo- 
dynamics and criticality of one-dimensional spin-1 bosons with strongly repulsive density-density 
and antiferromagnetic spin-exchange interactions. We analytically derive a high precision equation 
of state from which the Tomonaga-Luttinger liquid physics and quantum critical behavior of the 
system are computed. We obtain explicit forms for the scaling functions near the critical points 
yielding the dynamical exponent z — 2 and correlation length exponent v — 1/2 for the quantum 
phase transitions driven by either the chemical potential or the magnetic field. Consequently, we 
further demonstrate that quantum criticality of the system can be mapped out from the finite tem- 
perature density and magnetization profiles of the ID trapped gas. Our results provide the physical 
origin of quantum criticality in a ID many-body system beyond the Tomonaga-Luttinger liquid 
description. 

PACS numbers: 03.75.Ss, 03.75.Hh, 02.30.1k, 34.10. +x 



I. INTRODUCTION 

The study of spinor Bose gases is an active area of 
research in the field of cold atoms P, 0- In an opti- 
cal trap, the laser-atom interaction is determined by the 
induced electric dipole moment, thus the atoms are con- 
fined independently of their spin orientations. This has 
provided exciting opportunities of simulating quantum 
dynamics of spinor Bose-Einstein condensates in which 
the "vector" property of spinor atoms can be preserved. 
Several experimental groups have successfully demon- 
strated spinor BECs of 23 Na 0, i] and 87 Rb [fQ atoms 
in optical traps. In particular, the exquisite tunability 
with ultracold atoms confined to low dimensions has pro- 
vided unprecedented opportunities for testing the the- 
ory of one-dimensional (ID) exactly solvable many-body 
systems ■ These experimental developments have 

stimulated an extensive study of related exactly solvable 
models with <5-function interactions, see recent reviews 
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Quantum spinor gases with multi-spin states exhibit 
richer quantum effects than their single component coun- 
terparts. Spinor Bose gases with spin-independent short 
range interaction have a ferromagnetic ground state, i.e. 
the ground state is always fully polarized f2pL [2l| . In 
contrast to the two-component Fermi gases [22j | , the two- 
component spinor Bose gas with spin-independent s-wave 
scattering |2a - l25| has a ferromagnetic ground state as 
long as the interaction is fully spin independent. How- 
ever, ID spinor Bose gas with short-range density-density 
and spin-exchange interactions [H, [27j can display a dif- 



1.5 
1 

0.5 

X) 
CO 

i° 

-0.5 
-1 
-1.5, 



: \ M 

-\ \ 00 
\ \ v§.'4t'<W<ID 

: s\ 

-|| 


# spin +1 - 

• spin 
©spin-1 _ 


an©©© 


: V 



0.5 



H/8, 



1.5 



*xwcl05@physics. ami. edu.au 



FIG. 1: (Color online) Phase diagram in the [i-H plane show- 
ing the spin-singlet phase S of paired bosons, ferromagnetic 
phase F of spin-aligned bosons, and a mixed phase M of pairs 
and unpaired bosons. V stands for the vacuum. The S — M 
and F — M boundaries are determined by the critical fields (|8]) 
and 1)101) . respectively. The V — S and V — F boundaries are 
given in ([6]). The dashed-dotted lines are the extrapolation 
of the phase boundaries and ()12p in the strong coupling 
regime. 



ferent ground state, i.e., either a ferromagnetic or an 
antiferromagnetic ground state solely depending on the 
spin-exchange interaction. In this context, the spin-1 
spinor Bose gas with short-range delta-function interac- 
tion and antiferromagnetic spin-spin interaction is partic- 
ularly interesting due to the existence of various phases of 
quantum liquids associated with the Bethe Ansatz (BA) 
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solutions [26|, |28 30]. At zero temperature, this model 
exhibits three phases in the chemical potential - mag- 
netic field plane. These are (i) a spin-singlet phase of 
pairs of bosons with hyperfine states \F — l,mp = ±1) 
or two \F = l,mp = 0) bosons, (ii) a ferromagnetic 
phase of fully-polarized atoms in the hyperfine state 
\F = l,mp = 1) and (iii) a mixed phase of spin-singlet 
pairs and unpaired single atoms, see Figure [TJ 

Spinor Bose gases exhibit various phases of strongly 
correlated quantum liquids and are thus particularly 
valuable to investigate quantum magnetism and critical- 
ity. Near a quantum critical point, the many-body sys- 
tem is expected to show universal scaling behaviour in the 
thermodynamic quantities due to the collective nature of 
many-body effects. Thus a universal and scale-invariant 
description of the system is expected through the power- 
law scaling of thermodynamic properties [3l|, HH . Most 
recently, quantum criticality and universal scaling be- 
haviour have been experimentally investigated in low- 
dimensional cold atomic matter [33, 34]. These advances 
build on theoretical schemes for mapping out quantum 
criticality in cold atom systems |35l - l37l ]. In this frame- 
work, exactly solvable models of cold atoms, exhibiting 
quantum phase transitions, provide a rigorous way to 
treat quantum criticality in archetypical quantum many- 
body systems, such as the Gaudin-Yang Fermi gas [38| . 
the Lieb-Liniger Bose gas and a mixture of bosons 
and fermions [401 ] . 

Despite much work on the ID spin-1 bosons with 
strongly repulsive density-density and antiferromagnetic 
spin-exchange interactions [26|, [28|, [29[ , there has been no 
study of the quantum criticality of the model by using 
the exact solution. Contrary to the impression one might 
have, exact solvability does not guarantee that physical 
quantities of interest can be actually calculated by the 
BA solutions. The thermodynamic Bethe ansatz (TBA) 
equations for this model [29| involve an infinite number of 
coupled nonlinear integral equations that impose a num- 
ber of challenges to access the physics of the model. 

In the present paper, building on the method proposed 
in the stud y o f quantum criticality of the Gaudin-Yang 
Fermi gas [38| and the Lieb-Liniger Bose gas (3{|, we 
analytically study the quantum phase diagram, univer- 
sal thermodynamics and criticality of spin-1 bosons with 
strongly repulsive density-density and antiferromagnetic 
spin-exchange interactions. We derive a high precision 
equation of state of the system in experimental accessi- 
ble conditions, i.e., in the strong coupling regime and low 
temperatures. We also analytically derive the Tomonaga- 
Luttinger liquid (TLL) thermodynamics, quantum crit- 
ical exponents and universal scaling functions near the 
critical points associated with quantum phase transi- 
tions driven by the chemical potential and magnetic field. 
These scaling forms for the thermodynamic properties 
across the phase boundaries illustrate the physical origin 
of quantum criticality in this system, where the singular 
part of the thermodynamic properties involves a sudden 
change of density of state for either pairs or unpaired 



single atoms. 

The paper is organized as follows. In Sec[Hl we present 
the model and its corresponding TBA equations. In Sec 
IIII1 we analytically determine the phase diagram of the 
model at zero temperature. In Sec IIV1 we derive the 
equation of state and universal TLL thermodynamics in 
the physical regime, i.e., for strong coupling and low tem- 
peratures. In Sec [Vj we investigate quantum critical be- 
haviour driven by the chemical potential and magnetic 
field. The scaling functions near the critical points are 
obtained analytically. Sec|VT]is the conclusion. 



II. THE MODEL 

We consider N particles of mass m confined in ID 
to a length L with <5-intcracting type density-density 
and spin-exchange interactions between two atoms. The 
Hamiltonian is given by [l|, [26[ 



n = 



N 



i:) 2 



2—1 



c 2 Si ■ Sj]S(xi -Xj)+E z , (1) 



i<j 



where Si is the spin-1 operator with z-component (s 
1.0. I } . The interaction parameters cq = (go + 2g 2 ) /3 
and C2 = (.92 — <7o)/3 where gs = 4:irh 2 as jra. Here m is 
the particle mass and as represents the s-wave scattering 
length in the total spin S = 0, 2 channels. E z = —HS Z 
stands for the Zeeman energy, where H is the external 
field and S z the total spin in the z-component. In the 
above equation, we have set h — 2m = 1. 

Using the BA hypothesis, Cao et al. (2(| solved the 
model ((TJ) with antiferromagnetic spin-exchange interac- 
tion for c — Co = C2 > 0. The energy eigenspectrum is 
given in terms of the quasi- momenta {kj} of the parti- 
cles through E = Y^j= i tf i obeying the following set of 
coupled BA equations |26j 



N 



exp(ifcj-L) = TJ ' 



4c'i -r-r kj - A c 



2c'i 



kt - 4c'i 1 1 kj - A a + 2c'i 

a— 1 J 



N 

n 



A Q - k e + 2c'i 
A Q — kg — 2c'i 



M 

n 

0=1 



A, 



A^ + 2c'i 



A a - A 



p 



2c'i 



(2) 



Here d = c/4, j = l,...,N, a = 1,...,M and {A Q } are 
the rapidities for the internal spin degrees of freedom. 
The quantum number M is a conserved quantity obeying 
the relation M = N — S z . In this model, the antiferro- 
magnetic interaction leads to an effective attraction in 
the spin-singlet channel so that the singlet bosonic pairs 
comprise a spin singlet ground state. In the thermody- 
namic limit N,L — > oo, the sets of solutions {kj} and 
{A a } of the BA equations @ take a certain form, where 
the kj 's and A Q 's can form complex pairs kj = Xj ± id 
and Kj — Xj ± ic' where Xj is real. Notice that each 
pair of kj's share the same real part as a corresponding 
pair of Aj's. The bound states are associated with a pair 
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of \F = l,mp = ±1) bosons or two \F = l,mF = 0) 
bosons. In addition to that, we also have real kj's and A 
strings of the form A™ J = A"+i(n+l — 2j)c' , j = 1, . . . ,n 
describing spin wave bound states. 

At finite temperatures, the physical states become de- 
generate. The equilibrium state can be obtained by 
the condition of minimizing the Gibbs free energy G = 
E + E z — fiN — TS, where fj, is the chemical potential and 
S the entropy, see Yang and Yang's grand canonical de- 
scription [41] of the BA equations for the integrable Bose 
gas. The Zeeman energy Ez — —HS Z and entropy S are 
given in terms of the densities of charge bound states and 
spin-strings described above which are subject to the BA 
equations @. Minimizing the Gibbs free energy leads 
to a set of coupled non- linear integral equations, i.e., the 
TBA equations (see [29( for details) 



£l (k) = k 2 - n - H - Ta 4 * ln(l + e" 



-- 1 ' i . ' 



+T[ai - a 5 ] * ln(l + e 

OO 

-T^[a„_i + o n +i] *ln(l + e 



n=l 
„2 „/2 



--i(fc) 



£ 2 (fc) = 2(k J - c u - n) + T[ai - a 5 ] * ln(l + e ~ ) 

e 2 (fc) 

+T[a 2 - a 4 - a 6 ] * ln(l + e t ), 
0„(fc) = n + T[a„_i + a„+i] * ln(l + e t ) 



-T^T mn *ln(l + e-^). 



(3) 



n=l 



Here n = 1, 2, . . . , oo and the symbol * denotes the con- 
volution (/ * g{x)) — f(x — x')g(x')dx' , the func- 



tions a n 



and T nm are given in 



These 



7r (nc')' 2 -^?.' 2 

TBA equations are expressed in terms of the dressed en- 
ergies £i(fc), E2(k) and <f> n (k) for unpaired states, paired 
states and spin strings, respectively. They depend on 
the chemical potential /Lt, the external field H and spin 
fluctuations which are ferromagnetically coupled to the 
unpaired Fermi sea. 

The pressure per unit length of the system is given by 
p = P1+P2 with 



Pi 



T 
2^ 



ln(l 



-ei(k)/T 



P2 = ~ 
IT 



— 00 
oc 



ln(l 



-e 3 (fc)/T 



)dk, 



)dk, 



(4) 



corresponding to the pressures for unpaired bosons and 
spin-singlet pairs, respectively. 



III. PHASE DIAGRAM IN THE n - H PLANE 

The ground state properties and phase diagram at zero 
temperature can be determined by the dressed energy 
equations 



e 2 (k) 



f-i - H + a 4 * £i(k) + [a 5 - a{\ * e 2 (k), 



2(k 2 



„'2 



A 1 ) + [0-5 ~ ai] * £i(fc) 



+ [a 6 + a 4 - a 2 ] * e 2 (k) 



(5) 



which are obtained from the TBA equations ((3]) in the 
limit T — > 0. The negative part of the dressed energies 
e a {k), a = 1, 2 for k < Q a corresponds to occupied states, 
while the positive part of e a corresponds to unoccupied 
states. The integration boundaries Q a characterize the 
"Fermi surfaces" defined by e a (Q a ) — 0. In a canonical 
ensemble the bosonic pairs form a spin singlet ground 
state when the external field is less than a lower critical 
field. In this phase, the low energy physics can be char- 
acterized by a spin-charge separation theory of the U(l) 
TLL describing the charge sector and a Q(3) non-linear 
sigma model describing the spin sector [281 ] . However, 
if the external field exceeds an upper critical field, we 
have solely ferromagnetic single bosons with aligned spins 
along the external field. For an intermediate magnetic 
field, the spin-singlet pairs and spin-aligned bosons form 
a two-component TLL with magnetization [29J . However, 
in realistic experiments with cold atoms, ID systems can 
be realized by tightly confining the atomic cloud in two 
(radial) dimensions and weakly confining it along the ax- 
ial direction in a harmonic trap. Therefore, the phase 
diagram in the fi — H plane is essential for understanding 
quantum criticality of the trapped gas at finite tempera- 
tures. 

We may determine the phase boundaries by analysing 
the band fillings in the dressed energy equations (0 . The 
V — F phase boundary is established by the condition 
£%(k) < and £ 2 (fc) > 0. Then from equation ([5]) we 
have fid — —H. The V — S phase boundary is de- 
termined by £i(fc) > and £ 2 (fc) < that results in 
/i c2 = — eb/2, where = h 2 c 2 /(16m) is the binding en- 
ergy of the bound pair. For convenience we shall use the 
dimensionless units in the study of quantum criticality 
of the system, i.e., fl = /x/ej and h = H/e^. Thus the 
critical fields for the phase boundaries V — F and V — S 
read 



-h, 



A* c2 



(6) 



The (F — M) phase boundary is obtained by the re- 
quirement £i(±Qi) = and £ 2 (fc) < 0, yielding the set 
of equations 



1 f Ql Ei fx 1 ) 
ix{x) = 8x 2 -fi c3 -h+- / v dx', 

n J-Qi 1- + (x — x ) 



Mc3 h _ !_ f Ql 
8 8 8vr J_ Qi 1 + (Q x - x'Y 



h(x') 



dx 1 , (7) 



which give the critical field for the phase transition from 
a ferromagnetic phase of spin-aligned bosons into a mixed 
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phase of the pairs and unpaired bosons, 



1 2 



Mc3 



Qi 



£l(x) 



2 5tt 1 + 16a: 2 /25 
2 f® 1 ejjx) 

7T J_Q t 1 + 16x 2 



dx, 



rdx. 



(8) 



where Qi = Qi/c and £i(x) is given by 0. 

The S—M phase boundary is determined by the condi- 
tions £i(fc) < and £ 2 (±Q 2 ) = 0, from which we obtain 
the set of equations, 



e 2 (x) = 2 ( 8x z - - - nd 

2 I t J_Q a 



1 + (x - x') 



rdx 



Q-2 



i 2 (x') 



dx' 



3t J-q 2 1 + 4(x - x') 2 /9 

- 2 - r i2[xi) d^ 

IT J_q 2 1 + 4(x - x'Y ' 
1 Ac4 1 /" Q2 e 2 (a:') 



16 8 24tt J_q 2 l + 4(g 2 - x') 2 /9 



-dx' 



16tt 7_q 2 1 + (Q 2 - .x') 2 
1 Z^ 2 £ 2 (a;') 



8vr 7_q 2 l + 4(Q 2 -x') 



A2 



(9) 



that provide the critical fields for a phase transition from 
the spin-singlet phase of paired bosons into a mixed phase 
of pairs and unpaired bosons, 



for unpaired and paired bosons are presented explicitly 
in 

In the next section we will derive analytical expressions 
for the equation of state and universal TLL thermody- 
namics in the physical regime where t = T/e^ "C 1, i.e., 
for the strong coupling and low temperature regimes. 



IV. EQUATION OF STATE AND TLL 
THERMODYNAMICS 



The thermodynamics and the high precision of the 
equation of state of a system are the key informations 
that can be used to map out quantum critical phenom- 
ena and to make comparisons between theory and experi- 
ment. Recently, the equation of state of a two-c omp onent 
ultra-cold Fermi gas has been measured [42|, |43[ using 
theoretical schemes [44j. Such experimental advances 
provide exciting opportunities to test universal TLL and 
quantum critical phenomena in low dimensional many- 
body systems. 



if* s 2 (x) 



5tt J_q 2 1 + 16a; 2 /25 
J f® 2 £2 Or) , 

■rdX. 



t J-o, 1 + 16x 2 



(10) 



A. Equation of state 



Here Q 2 = Q 2 /c and £ 2 {x) is given by (HJ). 

In order to investigate quantum criticality of the sys- 
tem in the strong coupling regime, we need closed form 
expressions for the critical fields. By Taylor expansion of 
Eqs. (jH) and (fTO]) . we obtain the critical field values 



104 / 1 , 

t: \h--\ 11 

75tt 2 V ° ' 



1 8^2 / 1 
^ = ~2 + 15^^-2. 

, 32^2 [1 ,V 2912 fl \* ^ 



15tt V2 



225tt 2 \2 



which are in good agreement with the numerical results 
obtained from §8$ and (fTTJj) in the strong coupling regime. 
These asymptotic results (Eqs. flTT} and ([12])) can also 
be obtained by converting the critical fields obtained in 
the H — n plane [2^] into the \i — H plane, where the 
effective chemical potentials fjb\ = fi + H , fj, 2 = + £b/2 



The lack of analytic solutions of the TBA equations 
limits the ability to make physical predictions of the 
model at finite temperatures. In fact, the thermody- 
namic properties of the model at finite temperature are 
notoriously difficult to extract due to the presence of the 
bosonic nature and the spin-spin exchange interaction. 
Building on the method presented in (38| and consider- 
ing the physical region (strong coupling \c\ >• 1 and low 
temperatures) , we find that spin fluctuations are strongly 
suppressed by a strong field, i.e., H 3> T. Therefore we 
can analytically extract the spin wave bound state con- 
tributions to the unpaired dressed energy, see the third 
equation in ©. Moreover, we notice that the convolu- 
tion terms converge rapidly once £i, 2 (&) > in the TBA 
equations. Therefore, we are allowed to carry out a Tay- 
lor expansion with respect to c in the kernel functions a n 
in the TBA equations at low temperatures. Then, inte- 
grating by parts, we may obtain the dressed energies in 
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terms of polylogarithm functions up to order l/\c\ 
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FIG. 2: (Color online) Entropy vs temperature from TLL 
entropy (|30p and entropy calculated from the equation of state 
flgl-ITg]) for fi = -0.49 and /x = -0.45. The universal linear 
temperature dependent TLL entropy is broken down at the 
crossover temperatures t* which separate the TLL phase and 
the quantum critical regime, see Fig [3] 



Here Li s (z) — ^ZfcLi z k /k s is the polylogarithm function. 
The terms K = 8 Pl /(T\c\) and I n (z) = £~ =0 
are obtained from the so-called "string" or spin wave con- 
tributions. 

Using the above dressed energies and integrating by 
parts, we may calculate the pressure Q in a straightfor- 
ward way, with result 



Pi 



I>2 



Tifl 

2 

( 4lTh 2 \ 2 

I 2m/ 

2 

( 2ttH' 2 \ 3 
V 2m / 



pi 2m 3968^2 2m 

]c\*HF + i25|c| 3 "fiF 

3968pi 2m 181p 2 2m 



where we have denoted /. 
and 



7J 



2pi 

Id 



125|c| 3 S 2 108|c| 3 h 2 
' = Li s (~ e ^") w ith * 
16p 2 



(15) 

(16) 
1,2 



Ti 



A, 



2^¥ + 2/i 



/I 



5|c| 
1 

_2~7^ 

32pi 
5|c| 

' 3968 
1250F 



1984 



/I 



125V2tt 5 



(17) 



_P2_ 

3|c| 
/I- 



181 



ft 

108V2tt = 



(18) 



Thus an infinite number of TBA equations have been 
simplified to two coupled equations, making the thermo- 
dynamics of the model analytically accessible. 



The above expressions for the pressures (IT51) and fT| 
provide the precise equation of state of the system from 
which universal TLL thermodynamics and scaling func- 
tions near critical points can be further derived analyti- 
cally. We present a high precision equation of state in the 
later discussions of the singularities of thermodynamic 
properties near the quantum critical point as the tem- 
perature tends to zero. To evaluate physical properties 
we substitute Eqs. JT7J) and (JT5J) into Eq.([T5]| and ([Tol 
This provides two coupled equations for p\ and P2 , which 
can be solved by iteration. We discuss the quantum crit- 
icality of the system using these pressures in the next 
Section. 



B. Universal TLL thermodynamics 

In general, the free energy of a ID many-body sys- 
tem at low temperatures can be naturally attributed to 
low-lying excitations near the Fermi points. From calcu- 
lations of the finite temperature corrections to the free 
energy, one can extract universal TLL thermodynamics. 
This low energy phy sics can also be obtained from confor- 
mal field theory [45( . Here we further develop an efficient 
way to obtain the universal TLL thermodynamics in the 
mixed phase of the pairs and unpaired bosons from the 
TBA equations For temperatures T « 1 and the 
strong coupling regime |c| ^ 1, the gapless phase is in 
the region H ~ c 2 > T. Therefore, we can ignore the 
spin wave bound state contributions in this phase. The 
first few terms coming from an asymptotic expansion in 
the TBA equations are given in terms of l/|c| correc- 
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FIG. 3: (Color online) Contour plot of entropy S in the t — fj, 
plane from the equation of state (|15[) and (|16|l for two values 
of the external magnetic field (a) h — 0.49 and (b) h = 0.51. 
The dashed lines are the crossover temperatures determined 
in Figure [2] The intersection points between two dashed lines 
are the critical points in the /i — H plane, see Figure [T] (a) 
The dashed line separates the TLLs and TLLm from the 
quantum critical regimes, (b) The dashed line separates the 
TLLf and TLLm from the quantum critical regimes. The 
left-most-dashed line in (a) and (b) separates the vacuum V 
from the TLLs and TLLf, respectively. 



tions by 
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2pi 16p 2 
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5 c 3 c 



(19) 
(20) 



Higher order corrections can be calculated in a straight- 
forward manner (see the analysis in (|15p and (jlGI) ). How- 
ever, they are not necessary in the present discussion. 
From Eq. (jj) we obtain the pressures for spin-aligned 



single bosons and spin-singlet pairs 
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h 2 \~ 1/2 i "» 



e\de\ 



Pi = 



2m J it J l + e( £ ?-^?)/T : 
h 2 \- 1/2 V2 



2m 



where 



tt Jo l + e( £ §- A °)/ T ' 



2pi 16p 2 



A °"" + H + ^ Bid- 



A\ w 2n + 



h 2 c 2 32p 1 p 2 



2m 8 5|c| 



3|c! 



(21) 
(22) 

(23) 
(24) 



The integrals in (j2"Tj) and (|2"2"j) can be calculated explic- 
itly via the Sommerfeld expansion. We assume that there 
are two "Fermi seas", i.e., a Fermi sea of bound pairs 
with an effective chemical potential A\ and a Fermi sea 
of unpaired bosons with an effective chemical potential 
A\. From equations (|19H24[) and the relations n\ — J-^, 
n = n\ + 2n2 — Sj, by a cumbersome iteration we can 
obtain closed forms for the pressures 



Pi 



27r 2 n? 



+ 



P2 



+ 



3 

6ni 
3 

T 2 



1 



1 + 



3n 2 



6ni 

' w 

2ni 
48ni 
16n! 



96n 2 \ 
" 5|c| J 
32n 2 ' 

"2 



"2 

3]cl 



(25) 



(26) 



The Hclmholtz free energy per unit length is given by / = 
n/i — p. After a lengthy iteration, we obtain a universal 
leading temperature corrections to the free energy of the 
form 



f = fo- 



nT 2 



1 

vi 



1 

' ! 2 



(27) 



with /o the ground state already obtained in [29] . Here 

2(32n 2 - 10ni) 



and 



v\ = 2irni ( 1 



V2 = 7T712 



5\c\ 



2(48ni + 5n 2 ) 
15H 



(28) 



(29) 



are the charge velocities for unpaired and paired bosons. 
The entropy per unit length is given by s = — where 



7TT 

IT 



i 

i'2 



(30) 



We observe that in this gapless phase, spin wave bound 
state fluctuations are suppressed due to a strong exter- 
nal field. The suppression of spin fluctuations leads to a 
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universality class of a two-component TLL in the mixed 
phase of pairs and unpaired bosons, which we denote 
by TLLm- At low temperatures, the spin singlet phase 
persists as a single component TLL (denoted by TLLs) 
below a crossover temperature. The fully polarised sin- 
gle atoms can persist in a TLL phase (denoted by TLLf) 
below another crossover temperature. However, the TLL 
is not appropriate to describe quantum criticality, since 
it does not include proper thermal fluctuations for the 
quantum critical regime. In general the TLL persists be- 
low the crossover temperature at which the relation of the 
linear temperature-dependent entropy (or specific heat) 
breaks down. 

In Figure [5] we present the entropy as a function of the 
temperature using these two different approaches - the 
polylogarithm function result (I15M18[) and the Sommcr- 
feld expansion (f30f . The crossover temperature t* deter- 
mines the boundaries between the TLL regime and the 
quantum critical regime, see the contour plots of entropy 
in the t — fj, plane for two different values of the external 
field in Figure [3J The crossover boundaries are estab- 
lished by the points at which the TLL entropy (I3TJ1) breaks 
down, i.e., the entropy is no longer linear temperature- 
dependent. At finite temperatures, the system exhibits 
the characteristic V-shaped behaviour of quantum criti- 
cality. 



A. Criticality driven by a chemical potential 

Quantum critical behaviour is uniquely characterized 
by the critical exponents depending only on the dimen- 
sionality and the symmetry of the excitation spectrum. 
This is reflected by singularities in the thermodynamic 
quantities, such as density n, compressibility k = dnjd[i 
and magnetization M. They can be obtained from the 
derivatives of the pressure p with respect to \x and H. In 
order to identify universal scaling of the thermodynamic 
properties in the quantum critical regime, we will only 
take into account the first few terms in the equation of 
state (fT5j) and (fT6|) . To this end, the total pressure is 
simplified as 



(31) 



where we have defined a dimensionless pressure p 
p/\c\tb with the potentials 



-4i 



\i - 



2§x 



A 2 = 1 + 2/2- 



32pi 



16p 2 

5 ' 

_ P2 

3 ' 



(32) 
(33) 



and we have denoted the function /j = hi 

By iterating equations (I3T1) . (|3"2")l and 

the dimensionless density h = n/\c\, where 



-e T 
we obtain 



V. QUANTUM CRITICALITY 



4= f k l />Ai + V2f 1/2 A 2 ) , 



2^/2^ V 2 



(34) 



with 



Quantum criticality describes the critical behavior 
near a quantum phase transition, i.e., it describes col- 
lective behavior of a large number of interacting parti- 
cles at temperatures sufficiently low, such that quantum 
mechanics plays a crucial role in determining the distin- 
guishing characteristics |3lj . The quantum phase transi- 
tion occurs at absolute zero temperature as the param- 
eters of the system are varied. In the critical regime, a 
universal and scale-invariant description of the system is 
expected through the power-law scaling of thermodynam- 
ical properties. From the phase diagram of spin-1 bosons, 
see Figure [TJ we observe that the quantum phase transi- 
tion occurs as the driving parameters jl and h cross the 
phase boundaries at zero temperature. Although there is 
no true finite temperature quantum phase transition in 
a ID model, quantum criticality of a ID many-body sys- 
tem is associated with a universal crossover T* that sepa- 
rates the excitation spectrum from relativistic to nonrel- 
ativistic dispersion (4(||47j]. In the present spin-1 model 
of bosons, we can interpret the spin singlet phase as a 
TG gas of hard-core bosons with mass 2m and the spin- 
aligned ferromagnetic phase as a TG gas of single atoms 
with mass m. The mixed phase of two coupled TG gases 
is made of particles with mass m and 2m. 



a - i tl/2 f 1 + 1 ff 1 Y + 8tl/2 f* 

It / 2 \2 12* ! 2 

+ 15^V /l /2j + ^ 7f A/2/1/2. 

4+1/2 f . 2 +1/2 

1 + ^7=fbn - iz (fi ' 



A 2 = 



5V2tt 



t 



1 :,^ J] ~) T 120F^/ 2 



1447T 



fa)' 



JOlt 1 ,2 

75nV2 h/2h/2 ' 



The total density (fM)) depends on the density of single 
atoms and the density of paired atoms. Furthermore, 
using the standard thermodynamic relations, we obtain 
the magnetization M = M/\c\, 



M 




157T^ /ll/2/l2/2 y 



(35) 
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and the susceptibility x = \£b/\c\ 



m ^ 1/2 [l-^=f 1/2 + ^{f 1/2 



St 1 / 2 
'50F 



'If + *-(fi ? I 36t fi f 2 



1/2J1/2 



+ / - 1/2 ^l( /ll/2 ) }■ 



(36) 



By a lengthy calculation, the compressibility ft = «;e&/|c| 
is given by 



-f 1 



/2 



.1/2 



3*V2 i 6i i/ 2 
20^ /l/2+ 250F'' 1 / 2 



2V2tt" 
5^ 



/l/2(l 



2V2^ /l/2 60^ A / 2 



5V2tt 
120? 



i / Hi 1 / 2 ! 

^ 1 + T7T7^ / i/2 + 



10V27T 

/ 2 /2 1 + ^, 



25V27T ' 



/2 



2^2/! i/ 2 

40F 
40? 



(A2) 2 



^1/2 



/l/2 



.(37) 



These thermodynamic properties pave a way to ex- 
tract universal scaling functions in the vicinity of the 
critical points ft c . Following the procedure discussed in 
[38j , quantum criticality of these thermodynamic quanti- 
ties can be obtained in the limit T — > and T > \ft — ft c \ 
across each of the phase boundaries, with the results 



(V-F) 



n ~ — 
M ~ - 



=Lii C- 

7T 2 \ 



■e * 



4v / 2¥ 2 V / ' 



v/t T • / 2(fi-P c2 ) 

n ~ — 2^Lu ( — e * 
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h ~ — Ai0Lii ( 
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2(M-Mc3) 
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-e * 



(40) 



{S-M) 



n ~ n 04 - A 3 



0Lii ( 



— e * 



M ~ -A 4 0Lii f-e 



(41) 



Here M 3 = "03, n 04 , A«, with i = 1 ... 4, a and 6 are con- 
stants, independent of ft and t. They are given explicitly 



by 



with 



"-03 
"04 

Ai 
A 2 

A 3 
A 4 

a = ( h 



2vr02 V tt02 27r 2 y ' 

7T \ 67T 367T 2 / 



1 /'j 16^5 



20? V 5vrV^ 25tt 2 

4^ / 3^5 \ 
572vr 3 / 2 V 5ttV2/ ' 



1 



4V2tt 

1 

402^ 



320 8486 \ 
5tt + 75^ J ' 

160 86 \ 

57T 157T 2 / ' 



(42) 




In the above equations no3 and no4 are the background 
densities near the critical points /13 and u 4 , respectively. 
At quantum criticality, the above densities can be cast 
into a universal scaling form (see [U HH, [3(| ) , e.g., 



n(p,T) = n Q + T*- 



(43) 



Here the dimensionality d = 1 and the scaling function 
Cf(:r) = ALii(x) with a constant A. Consequently the 
dynamical critical exponent z — 2 and the correlation 
length exponent v = 1/2 can be read off from the univer- 
sal scaling form (|43|) . We observe that the spin-1 Bose 
gas belongs to the same universality class as spin-1/2 at- 
tractive fermions [HI due to the hard-core nature of the 
two coupled Tonks- Girardeau gases. 

In Figures [4] and [5] we plot the "scaled density" 
T~^ +1 ~~' {n{p, T) — no) versus ft for different values of 
the temperature near the critical points ft c i, ft c2 , ft C 3 and 
ftd- We observe that after an appropriate subtraction of 
the background density all curves at different tempera- 
tures intersect at the critical points, which is the hall- 
mark of criticality. In the regime of low polarization, i.e. 
P < P c , the true phase transitions from the vacuum into 
the spin-singlet paired phase and from the pure paired 
phase into the mixture of spin-singlet pairs and spin- 
aligned bosons occur as the chemical potential passes the 
lower critical point ft c2 © and the upper critical point 
ftd (fT2|) . respectively, see Figure [4] 

For large polarization, i.e., P > P c , the phase tran- 
sitions from vacuum into the ferromagnetic spin- aligned 
boson phase and from the spin-aligned boson phase into 
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FIG. 4: (Color online) "Scaled density" vs chemical potential 
for h = 0.49 at different temperatures T/e b = 0.005, 0.01, 
0.015, 0.02 and 0.025. The density curves at different tem- 
peratures intersect at the critical points. This feature can be 
used to map out the phase boundary V — S at jl C 2 = —0.5 
(Eq. in upper panel and map out the phase boundary 

S — M at /2 c4 ~ -0.489 (Eq. ([12])) in the lower panel. 
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FIG. 5: (Color online) "Scaled density" vs chemical potential 
for h = 0.51 at different temperatures T/et = 0.005, 0.01, 
0.015, 0.02 and 0.025. The density curves at different tem- 
peratures intersect at the critical points. This feature can be 
used to map out the phase boundary (V — F) at p, c \ — —h 
(Eq. in upper panel and map out the phase boundary 

(F - M) at /i c3 « -0.499 (Eq. |[TTJ|) in the lower panel. 



ibility across all phase boundaries, 



(V-F) k~ 

the mixture of spin-singlet pairs and spin-aligned bosons 

occur as the chemical potential varies across the lower ~ ^ 

critical point p, c i © and the upper critical point /x C 3 ~ 
(fTTj) , respectively, see Figure [5] The universal scaling be- 
haviour and the zero temperature phase diagram can be (F ~ -W) ^ — ^03 
identified from the finite temperature density profiles of 

the trapped gas where the local chemical potentials are (g _ k ~ K04 — — Li_i (— e ~ ^ ) . (44) 

replaced by the harmonic trapping potential. 




Here K03, K04 are the background compressibility in the 
Furthermore, we mention that similar calculations of vicinities of the critical points jl c ^ (Eq. (|lip ) and /i C 4 (Eq. 
the scaling function can be constructed for the compress- (|12l) ). whereas A5 and Xq are temperature- independent 
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constants 



K03 

A 5 



3\fa 3a 
4vrV2a V" ' tt\/2 
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(45) 
(46) 
(47) 
(48) 



Again, the critical exponents z — 2 and v = 1/2 can be 
read off the universal scaling function 3F(x) = ALi_i(x) 
in the universal form 



k{h,T) = k +T- +1 --F 



(49) 
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B. Criticality driven by a magnetic field 



0.02150 



0.48 0.50 0.52 0.54 0.56 0.58 0.60 

H/s h 

b 

FIG. 6: (Color online) Contour plot of the entropy S vs the 
external field H for fixed chemical potential p, = —0.495 in 
the T — H plane. The dashed lines are determined by com- 
paring the result from the equation of state (|15H16p and the 
TLL entropy (|30|l . The crossover temperatures separates the 
quantum critical regimes from the TLL phases. 



The quantum phase transitions driven by a magnetic 
field are particularly interesting. In the phase diagram 
Figure [TJ for fixed chemical potential we can vary the 
external field H to pass the phase boundaries (S — M) 
and (M — F). At finite temperatures, the three zero tem- 
perature quantum phases, i.e., the phase of singlet pairs, 
ferromagnetic phase of spin- aligned atoms and the mixed 
phase of pairs and single atoms, become the relativistic 
TLL of bound pair {TLL S ), TLL of single atoms (TLL F ) 
and a two-component TLL (TLLm) of paired and single 
atoms, respectively. We obtain the critical fields by con- 
verting the critical fields Eq. (H"2"j) and Eq. (jTTj) 



h c x = 
h c2 = 



32^2 
15tt 

1 /15tt 
4~ 



32 
45tt 2 



a+l) (50) 



' i+ 2 



These critical fields and the crossover temperatures 
can be observed in the contour plot of the entropy 
in the T — H plane, see Figure |51 The low energy 
TLL physics breaks down at the crossover temperature 
(dashed lines) where the dispersion of either bound pairs 
or unpaired single atoms becomes nonrelativistic. In par- 
ticular, in the vicinity of the quantum critical points h c \ 
and h C 2, the system exhibits two different quantum crit- 
ical regimes. From equations (IM1) . (|55|) and (|3"?| . we find 
the scaling function in the critical regime near the critical 
point JST 



where the constants are given by 
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In the vicinity of the quantum critical points h C 2, we 
obtain the scaling forms 



i ~ n 05 - A 7 \/tLii ( 
(S - M) <j M ~ -AgViLii (-t 
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(52) 
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where 
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with 



a 



157T 

5tt V~ 



/< - - j ■ (55) 

In this case, the background density Mqq — noe is equal 
to the total density at the critical point. The density (or 
magnetization) and compressibility can be recast into the 
universal scaling form 



M(h,T) = no + T^-W^!^), (56) 

k(j*,T) = K + Ti+ 1 -^^( ^ ih T ~J lc) y (57) 

with the same critical exponents as that for quantum crit- 
icality driven by the chemical potential, i.e. the dynam- 
ical critical exponent z = 2 and the correlation length 
exponent v = 1/2. 

In Figure [7] we show the magnetization as a function 
of the external field for different temperatures for a fixed 
chemical potential. All curves intersect at the critical 
point h c i without background magnetization. However, 
these curves intersect at the upper critical point h C 2 with 
appropriate subtraction of the background magnetiza- 
tion. We can obtain similar scaling behaviour for the 
densities near the critical points h c \ and h C 2 like that 
presented in Figures [4] and [5l It turns out that magne- 
tization can be used to map out the bulk phase diagram 
through the ID trapped gas at finite temperatures. 



VI. CONCLUSION 

Using the TBA equations, we have studied the quan- 
tum phase diagram, thermodynamics and quantum crit- 
ical behaviour of one-dimensional spin-1 bosons with 
strongly repulsive density-density and antiferromagnetic 
spin-exchange interactions. We have determined with 
high precision the equation of state from which the TLL 
thermodynamics, universal scaling functions and critical 
exponents have been obtained. The universality class of 
quantum criticality has also been discussed. 
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FIG. 7: (Color online) Magnetization vs external field for 
ji = —0.495. Upper panel: the intersection point of the 
magnetization curves at different temperatures gives the crit- 
ical external field /i c i at the boundary 5* — M. Lower panel: 
after a proper subtraction of the background magnetization 
Mo6 = n(t), the intersection point of the magnetization curves 
at different temperatures gives the critical external field h C 2 
at the boundary F — M. 



The scaling forms of the density, compressibility, mag- 
netization and susceptibility associated with the phase 
transitions driven by the chemical potential and magnetic 
field were rescaled to the universality class of quantum 
criticality of free fermions with critical exponent z = 2 
and correlation length exponent v = 1/2. It thus turns 
out that the quantum criticality of the spin-1 Bose gas 
belongs to the same universality class as spin-1/2 at- 
tractive fermions [38| due to the hard-core nature of the 
two coupled Tonks- Girardeau gases. We have also shown 
that the quantum criticality in ID systems involves a 
universal crossover from a TLL with linear dispersion to 
free fermions with a quadratic dispersion near the critical 
point. These scaling forms for the thermodynamic prop- 
erties across the phase boundaries illustrate the physical 
origin of quantum criticality in this system, where the 
singular part of the thermodynamic properties involves a 
sudden change of density of state for either pairs or un- 
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paired single atoms. The phase diagram, the TLL ther- 
modynamics and critical properties of the bulk system 
can be mapped out from the density and magnetization 
profiles of the trapped spinor gas at finite temperatures. 
Our results open the way to further study such universal 
features of ID many-body physics in experiments with 
ultracold atoms. 
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